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Abstract 

^ , In this paper we present the mathematical description and analysis ol a fractional-order regulated 

system in the state space. A little historical background of our results in the analysis and synthesis of 
the fractional-order dynamical regulated systems is given. The methods and results of simulations of 



the fractional-order system described by a state space equation equivalent to three-member fractional- 

s 



(^— N ' order differential equation with a fractional-order PD regulator are then presented. The possibility 

("■^ , of investigating the stability of such systems is also considered. 
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1. Introduction 



^ , Since the first references on fractional-order derivatives in the 17-th century, the theory of fractional-order 

! '~j ' derivatives and integrals was highly developed by many mathematicians. In the last five decades, many 

r^ [ authors made a great effort to apply this knowledge in practice. But only in the last decade can we 

j^ ' find some significant works concerned with the description, analysis, and synthesis of fractional-order 

regulated systems. 

Real objects are generally fractional-order, however, for many of them the fractionality is very low. 
But the control systems used so far were all considered as integer-order systems, regardless of the reality. 
Because of the higher complexity and the absence of adequate mathematical tools, fractional-order dy- 
namical systems were only treated marginally in the theory and practice of control systems, e.g. pi, 0, 0]. 
In these works the first generalizations of analysis methods for fractional-order control systems were made 
(s-plane, frequency response, etc.). 

At that time our research in this field began. On the basis of fundamental works [|l| Q we made 
at first an analysis of the fractional dimension of the controlled systems and properties of fractional 
regulators f^, §, 0, |§]- Until then we did not make difference between the terms fractal and fractional 
regulator. Very significant theoretical works appeared in 1993 and later |||, |lO[ ^ ^ |3| |lj|. This and 
the following results in this field were summarized in book Pq| . The fundamentals for practical utilization 
of numerical methods H, hO, ni|, based on the relation for approximation of the fractional derivatives 
Pq], and analytical methods |12{ for fractional-order systems simulation were made. But the proposed 
iterative numerical methods did not converge to the correct solution. The first explicit noniterative 
numerical methods for simulation of fractional-order systems were presented in works mTj Iq, M . In 




3|, |2|1 etc. After 
two exact methods 



works |lj, |l§l the first comparison was made of these noniterative numerical methods and corresponded 
analytical methods derived for fractional-order control system from analytical solutions published in |1J] . 
An example of experimental design of a fractional-order PD controller, with comparison of dynamic 
properties in fractional- and integer-order system with a fractional- and integer-order controller, designed 
for an integer-order system as the best approximation to given fractional-order system, was presented 
in works [|l^ |8| too. It followed from these results, that an application of integer-order regulator to 
fractional-order system is inadequate and with a change of system or regulator parameters can lead 
to system instability. This example was taken to many other works |l3, nsl 
experimental method of the fractional-order system parameters identification |1 

pi| , p2[ were derived. A great effort was devoted to elaboration of methods of fractional-order controllers 
synthesis. From purely experimental methods |17, 18 we continued with methods based on minimization 
of regulated square ||2^, gs) etc., and methods based on stability and damping measure |2j, g^, ^^ etc. In 
these works the fundamentals of fractional-order PI^D controllers synthesis are developed and now first 
works appear with application of these methods to the control of chaotic fractional-order systems p9, M . 
The development of the methods described in works J25, 24, 29,j29] required the elaboration of methods 
of fractional-order system analysis in the frequency domain |27|, |28| etc. The following Bode analysis was 
utilized for the fractional-order controllers synthesis and such regulated control systems stability analysis. 
This contribution is a continuation of the previous works and deals with the mathematical description 
and analysis of fractional-order regulated systems in the state space. We present results of simulations 
of a fractional-order system with the aim to investigate the stability of such system. 



2. Example of the fractional- order control system in state space 

For the definition of the control system we consider a simple unity feed-back control system illustrated 
in Fig.l, where Gs{s) denotes the transfer function of the controlled system and Gr{s) is the controller 
transfer function, both integer- or fractional-order. 
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Fig.l: Feed - back control loop 



The differential equation of the above closed regulation system for the transfer function of the con- 
trolled system Gs{s) = l/{a2s" + ais^ + ao) and the controller Gr{s) ~ K + TdS^ has the form 

a2 y^"^ (t) + ai j/(« (t) + Ta y^'^ (t) + (ao + K)y{t)^K w{t) + T^ w^''^ (t) (1) 

where a, f3, S, are generally real numbers and ao, ai, 02, K, Td are arbitrary constants. 

In works |31^, p3, E9| was presented a state space model that expresses the fractional-order derivatives 

x(")(t) = Ax(i)+B ii(t), 

y(t) ^ C x(t), t > 0. (2) 

This description is convenient for simple models of systems pl| with only one fractional-order derivation. 
In this contribution we propose state space model of the linear time invariant one dimensional system 
which expresses the first derivatives in the state space equations (0) and which has the classical state 
space interpretation for the fractional-order system too. On the right side of these equations we can then 
transfer more than one fractional-order derivatives of the state space variables (x^^^\t)). A disadvantage 
of this expression is that we cannot describe the state space equations in vector and matrix relations as 
in previous description (0). 

yit) = g{^^f^Ht),uit)), t>0. (3) 
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Fig. 2 : Unit step responses 



Fig. 3 : State trajectories 



We verified the above methods on an example from p7| ^q] . Assume the system described by differen- 
tial equation (|l|) with coefficients 02 = 0.8, ai = 0.5, aq = 1, a = 2.2, /3 = 0.9, X = 20.5, Td = 3.7343 and 
5 = 1.15. After its modification and with state space variables x{t) — xi{t),x'{t) — X2{t),x"{t) — X2{t) 
we can derive the following state space model equivalent to model O) 



x\{t)^X2{t) 



4it) = - 



ao+K (2_a) 



02 

y{t)^KxAt)+Tdx';-'>it) , t>0. 
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We can made other alternative state space models for the same system. The order of the integer-order 
system in the state space model was equivalent to the number of the state variables. In fractional-order sys- 
tems it does not hold. We propose to consider the order of the fractional-order systems according to the 
order of the highest derivative (integer- or fractional-order) in the resulting differential equation, e.g. (U^). 

Under condition x'(i) = we can obtain from state space model (0) the two equations of statics of this 
system, which are not algebraic equations as in integer-order linear systems, but differential equations, from 
which we can compute the coordinates of the equilibrium point, to which the state trajectories (H) tend. 

After discretisation of the first derivative by first differentiation in (H) we obtained the simple Euler 
methods for solving the state space model. For the fractional-order derivatives in (Q) we can take the 
relation from e.g. y, |l^, |l^ O, [l8|. In Fig. 2 and Fig. 3 is a comparison of the unit - step response of the 
classical numerical solution [[l7| , [iq] and the numerical solution of the state space model. The obtained 
state trajectories represent stable focal point for the above mentioned coefficients of the system (with 
stability measure St = —1.5 and damping measure Ti = 0.37, e.g. [^ ^) and unstable focal point for 
only one changed coefficient Td = 0.7343 of the same system. 

This work was partially supported by grant VEGA 1/7098/20 from the Slovak Agency for Science. 
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